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ABSTRACT 

We report the results of a CCD imaging survey, complimented by astrometric and spectro- 
scopic follow-up studies, that aims to probe the fate of heavy-weight intermediate mass stars 
by unearthing new, faint, white dwarf members of the rich, nearby, intermediate age open 
clusters NGC 3532 and NGC 2287. We identify a total of four white dwarfs with distances, 
proper motions and cooling times which can be reconciled with membership of these popula- 
tions. We find that WD J0643-203 in NGC 2287, with an estimated mass of M= 1 .02- 1 . 16M , 
is potentially the most massive white dwarf so far identified within an open cluster. Guided 
by the predictions of modem theoretical models of the late-stage evolution of heavy-weight 
intermediate mass stars, we conclude that there is a distinct possibility it has a core composed 
of O and Ne. We also determine that despite the cooling times of the three new white dwarfs in 
NGC 3532 and the previously known degenerate member NGC 3532-10 spanning ~90Myr, 
they all have remarkably similar masses (A/^0.9-lM Q ). This is fully consistent with the re- 
sults from our prior work on a heterogeneous sample of ^50 white dwarfs from 12 stellar 
populations, on the basis of which we argued that the stellar initial mass-final mass relation is 
less steep at Mi n i t >4M0 than in the adjacent lower initial mass regime. This change in the 
gradient of the relation could account for the secondary peak observed in the mass distribu- 
tion of the field white dwarf population and mitigate the need to invoke close binary evolution 
to explain its existence. Spectroscopic investigation of numerous additional candidate white 
dwarf members of NGC 3532 unearthed by a recent independent study would be useful to 
confirm (or otherwise) these conclusions. 

Key words: stars: white dwarfs; galaxy: open clusters and associations: NGC 2287; 
NGC 3532; NGC 25 16 



1 INTRODUCTION 

Only H, He and Li are believed to have been created in signifi- 
cant quantities during the earliest stages in the evolution of the 
Universe (e.g. IWagoner et al. Heavier elements, including 

the two main constitue nts of the Earth's atm osphere, O and N, are 
manufactured by stars (Burbidge et al.ll957b . These species are in- 
jected into the interstellar medium (ISM) of our Galaxy and oth- 



f E-mail:pdd@aao. gov.au 
© 2010 RAS 



ers primarily during the latter phases in the stellar lifecycle when 
stars, depending on their mass, evolve either towards a supernova 
explosion or up the asymptotic giant branch (AGB) and through 
the planetary nebula phase. Subsequen tly, they are mixed into the 
ambient gas dScalo & Elme green 2004) from which they will be 
incorporated into future generations of stars and planets. 

A detailed understanding of how stellar populations build up 
and modify the metallicity of the ISM over time is essential to fully 
comprehend the formation and the evolution of galaxies. An im- 
portant aspect of this is knowledge of the amount and the compo- 
sition of the gas a star of a particular initial mass (and metallic- 
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ity) will return to the ISM during it's lifetime. This is intimately 
linked to it's final evolutionary state. Most low and intermediate 
mass stars (A/i n i t ^>5-6MQ) are anticipated to end their lives as ei- 
ther He or CO white dwarfs while massive objects (Minit^lOM©) 
are expected to perish in Fe core-collapse Type II supemovae ex- 
plosions (SNe II). The fate of the heaviest intermediate mass stars 
(IMS) which lie in the interveni ng mass range, however, remains 
very uncertain (e.g. Siess 2006), despite these being comparable 
in number to the massive objects. Stellar evolutionary models sug- 
gest that the C in a partially degenerate core of an aged heavy- 
weight IM S ignites when it ac hieves M<;1.05Mq and burns to O 
and Ne dNomotol Il984l 1 19871) . before the star enters a phase of 
super-AGB evolution. The very high temperature environment at 
the bottom of the convective envelopes of super-AGB stars fa vours 
the p roduction of isotopes such as 7 Li, 14 N, 13 C and 26 Al dSiessI 
2010). If this chemically enriched envelope gas is removed before 
the core attains M~1.37Mq, the star likely ends lif e as an ultra- 
massive (M~l.l-l.35Mp)) O Ne white dwarf (e.g iRitossa et al.l 
ll99dlGarcia -Berroetal .Il997l) . However, if the degenerate interior 
reaches M~1.37M0 then electron captures onto 24 Mg and 20 Ne 
will lead to a core collapse and a weak typ e II electron capture su- 
pernova explosion (ECSNe; lNomotoll 19841) . The temperatures and 
particle densities inherent to ECSNe are anticipated to generate a 
somewhat more exotic cocktail of elements that also includes 64 Zn, 
70 Ge , 90 Zr and possibly heavy r-process elements (A>13Q; l0ianl 
l2008h. but which i s relatively poor in a-process elements and Fe 
dWanaio et all2009l) . 

Unfortunately, the treatment of key physical processes within 
these calculations remains rather rudimentary. This limits the ca- 
pacity of stellar evolutionary theory to make firm predictions re- 
garding the fate of stars in this mass range. For example, the 
details of how convective mixing is implemented in the models, 
especially during the He core burning phase, can systematically 
shift the initial mas s range of super-AGB stars by ~2Mq (e.g. 
iGil-Pons et ali r2007). It remains an open question as to whether 
during the super-AGB phase the atmospheres are cool enough to 
form dust (Poelarends et al. 2008). Dust appears to be essential for 
the efficient transfer of net outward momentum from the radiation 
field to the envelope gas , at least in thermally-pulsing AGB stars 
(e.g. lWachter et al.l2002l) . so it is still uncertain if single stars in this 
mass range can evolve into ultra-massive white dwarfs (UMWDs) 
at all. 

Investigations of field UMWDs have been unable to provide 
unambiguous evidence in support of an evolutionary link with 
heavy-weight IMSs, perhaps because this sizeable popula tion is 
heterogeneous in origin. For example, iLiebert et all d2005l) found 
that the spatial distribution of the highest mass white dwarfs de- 
tected in the extreme ultraviolet (EUV) sky surveys, including 10 
with M >1.1Mq, is remarkably similar to that of the young stellar 
populations in the vicinity of Gould's Belt. It has been argued that, 
given their space motions and cooling times, the UMWDs GD 50 
and PG0136 +251 could be ass ociated with the young Pleiades 
open cluster dDobbie et al.| [2006). In contrast, however, the empir- 
ical white dwarf mass distribution, when compared to the output 
of traditional population synthesis calculations, suggests that the 
majority of UMWDs may be produced through the merging of the 
lower mass components of pri mordial close binary systems, at least 
at the higher galactic latitudes dYuanll 19921 ; ILiebert et al.l2005h . 

The study of white dwarfs that are demonstrably bound mem- 
bers of young open clusters can provide better insight into the 
fate of heavy-weight IMSs. This is because these populations are 
co-eval and have comparatively well determined ages, making 



it possible to place tighter constraints on the prior evolution of 
their degenerate members. With ready access to large telescopes, 
the number of high mass white dwarfs known to be members 
of open clusters has increased substantially during the last few 
years (~30 degenera tes wi th M>0.8Mr in ~10 populations, e.g . 
IWilliams et alf2004l2009l : lKalirai et al]2005l ; lDobbie et alj2004r) . 
However, there remain very few, if any, known cluster white dwarfs 
with M<;1.05Mq. Amongst the most pr omising candidates are 
LB 1497 (Pleiades ; iBergeron et all 1 19921) . KRR 24 (NGC 2099; 
i Kalirai et alj|2005t) and four objects in NGC 2168 dWilliams et ail 
2009), but their estimated masses are either borderline or have large 
uncertainties due to the rather modest quality of the spectroscopy 
used to determine their fundamental parameters. It is unclear if this 
lack of compelling candidate cluster UMWDs is a manifestation 
of a real deficit or merely down to counting statistics. At present, 
NGC 2168 is by far the richest young population in which the white 
dwarf cooling sequence has been spectroscopically studied down 
to its termination. Extrapolation of the cluster mass function be- 
yond the main sequence turn-off mass suggests that the central 
regions of NGC 21 6 8 should conta i n ~18 white dwarf members 
dKalirai et al]|2003h . lWilliams et aid d2009l) identified 12-14 possi- 
ble degenerate members from their UBV imaging of the cluster 
core which is marginally consistent with this prediction (P~0.09 
for 12 objects), especially considering that a non-zero proportion 
of the cluster white dwarfs will be hidden fr om view by bin ary 
companions in earlier evolutionary phases (e.g. lWilliamsl2004h . 

In light of this we have undertaken a new survey of 
two open clusters, NGC 3532 and NGC 2287, that aims 
to identify and characterise the white dwarf remnants of 
their recently deceased heavy-weight IMS members. Both 
targets are rich, nearby, young/intermediate age popula- 
tions (r=300±25Myr and T=243±40Myr, respectively , 
Mevnet et al.lll993l; iKoester & Reimersl 1 19931 : lHarris et al.l 1 19931 : 
Kharchenko et al J 120051 : ISharma et alj 120061) which, as outlined 
in dDobbie et al.l l2009h . have been previously searched photo- 
graphically for white d warfs (e.g. Romanishi n & Angell Il980t 
iReimers & Koesteilll989l) . However, none of the six photometric 
candidates subsequently confirmed through sp ectroscopy as white 
dwarf mem b ers of these two populations dKoester & Reimersl 
Il98ll 1 19931 ; iDobbie et al] l2009h appear to have formed from 
progenitors with M^5Mq, despite both clusters being relatively 
rich. It is possible the remnants of these stars have remained 
undetected until now because the y lie beyond the d etection limits 
of the photographic plates (e.g. Idem et all 1201 ll) . Indeed, the 
total mass of NGC 3532 is estimated to be M to t >2000M Q 
dFernandez & Salgadolfl980l) . This is approximately t wice that of 
the 150Myr old cluster, NGC 2516 dJeffries et al.l200ll) . the central 
regions of which harbour fo ur massive white dwarfs d escended 
from stars with M^5M dKoester & Reimersl 1 19961) . On the 
assumption that the initial mass functi ons of these populations 
were of a similar power law form (e.g. ISalpeterlll955l) . it could 
be expected that 4±2 degenerate remnants from stars with 
A/<;5.5Mq reside within NGC 3532 and approximately half this 
number within NGC 2287. 

In subsequent sections we outline our new, more sensitive 
charged coupled device (CCD) based survey of these clusters and 
detail our initial selection of candidate white dwarf members. We 
describe our spectroscopic follow-up observations of these sources, 
the analysis of the spectral datasets and the determination of the 
fundamental parameters of the new white dwarfs uncovered. In 
Section [4] we estimate the distances and proper motions of the de- 
generates and compare these to the properties of their proposed host 
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Table 1. Details of the four white dwarf members of NGC 2516 derived from FORS1 (top) and FORS2 (bottom) observations. The tabulated 
spectroscopic signal-to-noise estimates correspond to per resolution element over the range A=4150-4300A. Masses and cooling times for each 
star have been estimated using the mixed CO core composition "thick H-layer" evolutionary calculations of the Montreal Group (e.g. Fontaine, 
Brassard & Bergeron 2001). The quoted limits on the near-IR magnitudes of the white dwarfs correspond to 3<r detections in the images. 
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Figure 1. The results of our fitting of synthetic profiles (thin black lines) to 
the observed Balmer lines, H-/3 to H-8 , of the four white dwarf candidate 
members of NGC 25 16 (thick grey lines). The flux^ units are arbitrary. 



clusters. For those white dwarfs we find to be probable members of 
NGC 3532 and NGC 2287, we derive progenitor masses and con- 
clude by discussing them in the context of stellar evolution and the 
stellar initial mass-final mass relation. First, however, we briefly re- 
visit the four high mass white dwarf members of NGC 25 16, which 
have traditionally been used to define the upper end of the IFMR. 
We use our new and improved spectroscopy to place tighter limits 
on their fundamential parameters and inferred progenitor masses. 



2 THE NGC 2516 WHITE DWARFS REVISITED 

The open cluster NGC 25 16 was recognised several decades ago 
as an attractive target f or investigating the form of the upper 
end of the IFMR (e.j 



(e.g. 140Myr, Mevnet et al 
120Myrs, iKharchenko et al 



Reimers & KoesteJ Il982h. It is young 
Sung et alj I2OO2E 



158Myrs, 

compa ratively nearby, (e.g . 
430 + ggpc. Tpachs & Kabudl 19891 : 407+ggpc. lTerndrup et al.ll2002t 
404lg pc. I An et alj 12007) and has a met allicity close to solar 
(e.g. [Fe/H]=-0. 10± 0.04, Isung et al J I2OO2I ; [Fe/H]=+0.01±0.07, 
iTerndrup et al.lE002l) . Moreover, with a mass that is approximately 
twice that of the Pl eiades it is relatively rich (M to t=1240-1560M Q , 
I Jeffries et al.l200ft) . A deep UV and red photographic survey of 25 
square degrees of sky cen tered on the cluster was performed by 
iReimers & Koester] jl982h and led to the identification of nine blue 
candidate white dwarf members. Ultimately, four of these were 
confirmed through spectroscopy to be DA cluster members and 
found to have masses significantly greater than the objects populat- 
ing the peak of the field white dwarf mass dist r ibution (NGC 25 16- 
1, -2 , -3 and -5, IReimers & Koested 1 19821 ; I Koester & Reimersl 
1996). Although early optical spectra, obtained with the European 
Southern Observatory's (ESO) 3.6m telescope, were of sufficient 
quality to allow confirmation of the trend that higher mass stars 
evolve to become more massive degenerates and to argue that the 
maximum mass of a whi te dwarf progenitor is Mi n i t <;6MQ (e.g. 
iKoester & Reimeri 1 1996b , their rather modest signal-to-noise at 
A5s4000A, suggests significant further reduction of the statistical 
uncertainties on the mass and cooling time estimates for these ob- 
jects could be achieved with modern data. 

Therefore, we re-observed the white dwarf members of 
NGC 25 16 with ESO's 8m Very Large Telescope (VLT) UT2 and 
the FOcal Reducer and low dispersion Spectrograph (FORS1). A 
detailed description of the FORS instruments may be found on the 
ESO webpageo The data were acquired in service mode within the 
periods 2007/04/24-27 and 2007/10/06-11/21. Since these targets 
are comparatively bright, we specified fairly relaxed constraints 
on the sky conditions and thus these observations were generally 
undertaken in poorer seeing and/or with some cloud present. All 
the FORS1 data were acquired using the 2x2 binning mode of 
the E2V CCD, the 600B+12 grism and a 1.6" slit which gives 
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Table 2. Inferred progenitor masses of all the known white dwarf mem- 
bers of NGC 2516 for the most likely cluster age (column 2) and plau- 
sible limiting ages (columns 3 and 4). The white dwarf cooling times 
were derived using CO core tracks and the stellar lifetimes have been 
tr ansformed to masses using the solar metallicity evolutionary models 
of lGirardi et"aH feOOOl) . 
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Figure 2. A l'x 1' J band image centered on NGC 2516-2 (N top, E left). 
A weak source is seen within ~1" of the position of this wh ite dwarf as 
measured from photographic plates (Reim ers & Koesterll982l) . 



a notional resolution of A/AA~500. Flat and arc exposures were 
obtained within a few hours of the acquisition of each of the sci- 
ence frames. The spectra were reduced and extracted following the 
procedure outlined in Section 13.31 A spectrum of the featureless 
DC white dwarf WD0000+345 was obtained with an identical set- 
up during this period and was used to remove the remaining in- 
strumental response. We have used a grid of synthetic spectra and 
followed the procedures described in Sections |4. H and [4.21 to re- 
determine the effective temperatures and surface gravities of these 
four stars (Table [T). To serve as a check on the repeatability of the 
output from our spectral analysis procedure, NGC 2516-1 was re- 
observed by us with the VLT UT1 and FORS2 on 2010/02/06 (see 
Section [3~3t , The excellent agreement between the parameter esti- 
mates for NGC 2516-1 we derive from the two distinct sets of ob- 
servations, lends confidence to our reduction and analysis process. 
Subsequently, for internal and external consistency, we have used 
the CO core, thick H-layer evolutionary models of iFontaine et al.l 
l l200ll) to determine their masses and th eir cooling times (Ta ble[T|>, 
From the earlier ESO observations, Ko ester & Reimersl (fl996) 
determined that NGC 2516-2 and NGC 2516-5 most likely have 
MM.OMp), with a best estimate for the former of M=1.05M©. 
IClaveretal.ld200ll) used these results and the stellar evolutionary 
models of dGirardi et al.ll2000l) to infer the masses of their progen- 
itor stars to be M; n i t ~6M0. Since these latter estimates are con- 
sistent with the initial masses at and above which stars of solar 
metallicity are predicted by that generation of the Padova group 
calculations to reach sufficiently high central temperatures to ig- 
nite C and, also, the former values are larger than the H-depleted 
core masses in these models (A/=O.97M0), it has been suggested 
that the N GC 25 16 wh ite dwarfs may be composed of O and Ne 
(e.g. IWeidema^n1l2005h . However, based on our new, improved, 
spectroscopy we find that it is most likely that all four (known) 
white dwarf members of the cluster have M<1.0Mq. This places 
them below more recent theoretical predictions for th e minimum 
mass of ONe core wh ite dwarfs (e.g. M~ 1 .05- 1 . 1M , Isiessl2007l ; 
iGil-Pons et alj|2003h . All four of these objects also now appear to 
lie within a fairly narrow (final) mass range (Table [T) despite ex- 
tending over ~30Myr in the cooling age domain (corresponding to 
~0.5Mq variation in initial mass e.g. Table [2}- Some stellar evo- 
lutionary models predict that the masses of ONe remnan ts are a 
strong function of initial mass (e.g. Ilben & Tutukov|[l985j) . More- 
over, at the most probable age for this cluster, the inferred progeni- 



tor masses are now all some way below 6Mq . So, in light of these 
new observational results and more recent theoretical predictions 
we propose that the cores of the NGC 2516 white dwarfs are com- 
posed of C and O rather than O and Ne. If there are ONe core 
UMWDs within this cluster, they remain to be found. 

As the presence of a close companion can substantially impact 
the evolution of a star, it is also worth briefly discussing here the 
new near-IR imaging (J and Ks) of these objects recently obtained 
as part of the ES O Visible and Infrared Survey Telescope for As- 
tronomy (VISTA, lEmerson et ai1l2006l; IPalton et al.l l2006) Hemi- 
sphere Survey. This reveals no evidence of unresolved cool com- 
panions with A / />0.04M Q to NGC 2516-1, -3 or -5. However, this 
data hints at a possible near-IR excess to NGC 25 16-2 (see Figure|2] 
and Table [T). These observations were obtained during Decem- 
ber 2009 in modest seeing (~l"). Assuming the source detected 
at 07 h 57 m 50.79 s -60°49'55.1" (J2000.0) in the J band imaging 
is associated with NGC 2516-2 (there is no detection in the shal- 
lower Ks imaging, see Table QJ, then we appear to be observing 
about a factor two too much flux in this waveband (we estimate 
a < 10% likelihood of a chance alignment based on the density of 
sources of similar or greater brightne ss in the image). According 
to the grids of synthetic photometry of lHolberg & Bergeronl (2006) 
a T e ff =35000K, log ,g=8.5 white dwarf could be expected to have 
V-J=-Q.l\ (M./=11.2) but from the VISTA photometry we esti- 
mate that NGC 2516-2 has V- J=0.19±0.20. Although we have not 
applied any transformations to either the synthetic or the observed 
photometry to account for the slight differences in the filter sys- 
tems, any discrepancies are expected to be very small. Moreover, 
foreground reddening towards th e cluster is relatively low , YL{V- 
J)^0. 16, based on E(B-V)=0 .07 jKharchenko et alj2005h and the 
interstellar extinction curve of lFitzpatrickl(ll999l). We note that ac- 
cording to the NEXTGEN models of iBaraffe et all d 19981) . an ob- 
ject with M./~ll (Mk s ~10) has a mass of only ~50Mj up itor at 
the age of NGC 25 1 6. Deeper J and K data for NGC 25 1 6-2 should 
be obtained, in good seeing, to confirm or otherwise the presence 
of a substellar companion. 
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Figure 3. Schematic plots of the open clusters NGC 3532 (left) and NGC 2287 (right) with the regions in each which have been surveyed with the ESO 2.2m + 
WFI, highlighted (large open rectangles). Candidate white dwarf members of the clusters identified from their location in B—V,V colour-magnitude diagrams 
(five-point stars, both solid and open) and objects which have been followed-up spectroscopically with the VLT + FORS (open circles) are overplotted. The 
locations of previously known white dwarf members of the two clusters are also shown (squares). 



3 OBSERVATIONS AND SELECTION OF TARGETS 
3.1 The optical imaging surveys 

We retrieved B and V band imaging covering ~0.3 square de- 
grees of sky towards each of NGC 2287 and NGC 3532 from the 
ESO data archive^ see Figure [5}. These observations, which con- 
sist of a short (20 or 30s) and two 240s integrations per filter, were 
undertaken as part of the ESO Imaging Survey on the nights of 
1999/12/03 (dark) and 20 00/02/24 (bright) r espectively, with the 
Wide Field Imager (WFI; iBaade et al.ll l999h and the 2.2m tele- 
scope, located at La Silla, Chile. The ESO WFI consists of a mo- 
saic of eight ££1/4096x2048 pixel CCDs and covers an area of 
34' x33' per pointing. The data were reduced using the Cam bridge 
Astro nomical Survey Unit CCD reduction toolkit dlrwin & Lewis! 
l200ll) . We followed standard procedures, namely, subtraction of the 
bias, flat-fielding, astrometric calibration and co-addition. Subse- 
quently, we performed aperture photometry on the reduced images 
using a circular window with a diameter of 1.5 x the full width half 
maximum of the mean point spread function (~ 1 .0- 1 .3")- Finally, 
we morphologically classified all sources detected in the combined 
frames and constructed merged catalogues of photometry. 

As these archival datasets are more than a decade old, detailed 
sky transparency information (e.g. LOSSAM) for the time at which 
these observations were conducted is not available to us. Accord- 
ingly, to gauge the sky conditions we have examined the instru- 
mental magnitudes of a series of bright (but not saturated) stars in 
the three sequential sub-integrations per filter. These were found 
to exhibit only small systematic differences of the order of a few 
percent from frame to frame. We have taken this to indicate that 
the sky was largely clear at the time of these observations. Subse- 
quently, to convert the B and V instrumental magnitudes onto the 

2 http://archive.eso.org/cms/ 



Table 3. A summary of the distances and reddenings we adopted in our 
photometric search for candidate white dwarf members of NGC2287 
and NGC3532 and the cluster ages and metallicities we assumed in es- 
timating progenitor masses. 



Cluster 


V-M v 


E(B-V) 


t (Myr) 


[Fe/H] 


NGC 2287 


9.30 


0.01 


243±40 


0.00 


NGC 3532 


8.35 


0.05 


300±25 


0.00 



standard Johnson system, we adopted zeropoints, colour terms and 
airmass co-efficients of k =24.53 and 24.lfl ki=0.31 and -0.09 
and k2=0.20 and 0. 1 10, respectively. While we have not performed 
a rigorous photometric calibration of these data, we believe the ac- 
curacy of our magnitudes, where sufficient photons have been de- 
tected, is good to ~10%. This is sufficient for the purpose of our 
initial surveys. We estimate from the number of point-like sources 
detected in the combined frames that these datsets are <;50% com- 
plete at 73,1/^21.5 and B, V?s22 for NGC 3532 and NGC 2287 re- 
spectively. 

3.2 Selection of candidate white dwarf members of 
NGC 3532 and NGC 2287 

We have used our catalogues of photometry to construct B-V,V 
colour magnitude diagrams for all point-like sources in the opti- 
cal data (Figure 0- As a guide to the location of the white dwarf 
cooling sequence of each cluster we have highlighted the previ- 
ously known degenerate members within the surveyed areas (filled 

3 www.eso.org/sci/facilities/lasilla/instruments/wfi/inst/zeropoints/ColorEquations/ 

4 www.eso.org/sci/facilities/lasilla/telescopes/dlp5/misc/Extinction.html 
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Figure 4. B — V,V colour-magnitude diagrams for all objects classified as stellar (dots) in the ESO WFI imaging of NGC3532 (left) and NGC2287 (right). As 
a guide to the likely location of the cluster white dwarfs we have highlighted the pr eviously known clust er white dwarfs recovered here (squares) and have also 
overp lotted evolutionary tracks for O.7M0, I.IMq (both CO core, thick H-layer; iFontaine et al .1200 lh and 1.1 6Mq (ONe core, thick H-layer; Althau set alj 
2007) white dwarfs. These tracks have been shifted to distances of V-My=835 and V-Mv=930 and adjusted to allow for a foreground reddenings of 
E(B-V)=0.05 and E(B-V)=0.01 for NGC 3532 and NGC 2287 respectively. The remaining symbols to the objects are as defined in the caption of Figure|3] 



squares). Additionally, we have overplotted evolutiona ry tracks for 
M=0 .7M W , M=1.1M (CO core, thick H-laver: IFontaine et alJ 
2001 ) and M=1.16M Q (ONe core, thick H-laver; lAlthaus et all 
20071) white dwarfs. These have been modified to account for the 
distance to and foreground reddening towards each cluster i.e. V- 
M v =8.35, E(B-V)=0.05 for NGC 3532 and V-M v =9.30, E(B- 
10=0.01 for NGC 2287 (see Table [Hi. Re-assuringly, a glance at 
Figure [4] reveals that the previously known white dwarf cluster 
members recovered here lie close to these theoretical tracks. 

Since we are primarily interested in identifying the oldest, 
most massive cluster white dwarfs, for each population we have 
used the V magnitude of the faintest previously known degenerate 
member for NGC 3532 and V«20.5 for NGC 2287) as 

a guide as to where to set the bright limits of our candidate selec- 
tion boxes (dashed lines). Moreover, as the lower end of the theo- 
retical tracks plotted in Figure |4] correspond to the expected loca- 
tions of white dwarfs which have cooled for a time approximately 
equal to the maximum likely age of each cluster (r~330Myr for 
NGC 3532 and r«280Myr for NGC 2287), we have used these to 
define the faint limits of our selection boxes. In determining our 
selection criteria with respect to colour, we have adopted as the 
red limit the coolest point of the three evolutionary tracks plotted 
in Figure [4] This represents a M=0.TMq white dwarf at the max- 
imum probable age of each cluster. The blue limit is defined by 
the colour of the A/=1.16Mq track and includes an allowance for 
the estimated uncertainty in the B-V colours (~0.15 mag.). The 
application of these selection criteria leads to the identification of 
seven and four candidate white dwarf members of NGC 3532 and 
NGC 2287, respectively. Summary details for these objects, includ- 
ing their co-ordinates and provisional magnitudes, are given in Ta- 



ble E] We note that NGC 2287-WDC J0646-2032 has been previ- 
ously identified as a potenti al white dwarf member of NGC 2287 
dRomanishin & Angei|[l980l) . 



3.3 Follow-up spectroscopy with the VLT and FORS. 

We have obtained follow-up low resolution optical spectroscopy 
spanning A?s3750-5200A for six of the candidates listed in Table|4] 
using the ESO VLT (UT1) and FORS2. These data were acquired 
in visitor mode on the nights of 2010/02/06-07. While seeing was 
generally good throughout this period, sky transparency on the first 
night was less than optimal with long spells of moderately thick cir- 
rus. Fortunately, during the second night the sky was largely clear. 
The data were captured using the 2x2 binning mode of the blue 
optimised E2V CCD, the 600B+24 grism and a 1.3" slit, which 
provides a notional spectral resolution of A/AA~600. Flat and arc 
exposures were obtained within a few hours of the acquisition of 
the science frames. 

The CCD data were debiased and flat fielded using the I RAF 
procedure CCD P RO C. Cosmic ray hits w ere removed using the rou- 
tine LACOS SPEC ( Ivan Dokkuml200lh . Subsequently the datasets 
were extracted using the APEXTRACT package and wavelength 
calibrated by comparison with He+HgCd arc spectra. Remaining 
instrumental signature was removed using a spectrum of the fea- 
tureless DC white dwarf LHS 2333 which was obtained using an 
identical set-up towards the beginning of our programme. The ap- 
proximate S/N of the co-added datasets is listed in Table|5] 

In error, we also spectroscopically observed a source in the 
direction of NGC 2287 which our photometry indicates is not asso- 
ciated with this cluster (concentric open circles in Figures [3] right 
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Table 4. A summary of the astrometric and photometric properties of the eleven candidate white dwarf members of NGC 3532 and 
NGC 2287 unearthed by our ZJ, Vsurvey imaging of 0.3 sq. degrees of sky towards each cluster. In column 2 we include the ID number 
of candidates towards NGC 3532 as listed in Table 2 of lClem et all <201 II) . 



Survey 


Other 


RA 


DEC 


V 


B-V 


VLT+FORS 


Designation 


IDs 


J2000.0 






? 






m s \ 


(°"\) 








NGC 3532- WDC Jl 105-5857 


WD Jl 105-585, 198279 


11 05 23.86 


-58 57 23.00 


20.6 


0.01 


Y 


NGC 3532- WDC Jl 106-5856 


161574 


11 06 16.87 


-58 56 04.90 


20.0 


0.15 


N 


NGC 3532- WDC Jl 106-5905 


WD Jl 106-590, 160608 


11 06 18.35 


-59 05 17.33 


20.1 


0.07 


Y 


NGC 3532- WDC Jl 106-5843 


WD Jl 106-584, 134356 


11 06 51.65 


-58 43 49.03 


20.2 


0.11 


Y 


NGC 3532- WDC Jl 106-5847 


128702 


11 06 59.07 


-58 47 33.85 


20.8 


-0.03 


N 


NGC 3532- WDC Jl 107-5842 


WD Jl 107-584, 109695 


11 07 21.97 


-58 42 12.54 


20.5 


-0.36 


Y 


NGC 3532-WDC Jl 107-5848 


104388 


11 07 28.74 


-58 48 24.25 


20.8 


0.07 


N 


NGC 2287-WDC J0644-2033 




06 44 39.84 


-20 33 34.08 


21.0 


-0.25 


N 


NGC 2287-WDC J0645-2025 


WD J0645-202 


06 45 22.07 


-20 25 09.67 


21.0 


-0.12 


Y 


NGC 2287-WDC J0645-2050 




06 45 29.93 


-20 50 18.45 


21.0 


-0.32 


N 


NGC 2287-WDC J0646-2032 


WD J0646-203, NGC 2287-4 


06 46 01.08 


-20 32 03.42 


20.9 


-0.27 


Y 




J1 106-584 




J1 107-584 




T eff =20923+™=K log g=8.61+££ 



-100 

AX (A) 



AA (A) 



J0645-202 




J0646-203 




-100 

AA (A) 





AA (A) 



Figure 5. The results of our fitting of synthetic profiles (thin black lines) to 
the observed Balmer lines, H-/3 to H-8 , of the six white dwarf candidate 
members of NGC 3532 and NGC 2287 (thick grey lines). The flux^ units 
are arbitrary. 



and[4] right). We include this object in the remainder of our inves- 
tigation as a test of our other membership selection criteria. 



4 FUNDAMENTAL PARAMETERS OF THE WHITE 
DWARF CANDIDATE MEMBERS OF NGC 3532 AND 
NGC 2287 

The energy distributions of all spectroscopically confirmed white 
dwarfs display the strongly pressure broadened H-Balmer lines 
characteristic of the spectral type DA. To obtain estimates of the 
fundamental parameters of our six candidates (and our test object) 
so that we can further scrutinise their cluster membership status, 
we have compared the Balmer lines in each spectrum to synthetic 
line profiles obtained from model atmosphere calculations. 



4.1 The model atmosphere calculations 

The grid of synthetic pure-H spectra we use here is virtually iden- 
tical to the one used in our previous work on these open clus- 

I ■ II ~TJ 

ters (i.e. Dobbie et al. 2009) with the exception that it has been 
extended to slightly lower effective temperatures (T e jj=12500- 
35000K). To re-iterate, our models have been generated using 
relatively recent versions of the plane-parallel, hydrostatic, non- 
local thermodynamic equilibrium (non-LTE) atmosphere and spec- 
tral synthesis codes TLUSTY (v200; Hubeny 1988, Hubeny & 
Lanz 1995) and SYNSPEC (v49; Hubeny, I. and Lanz, T. 2001, 
http://nova.astro.umd.edu/!. The model H atom employed in the 
calculations incorporates the 8 lowest energy levels and one su- 
perlevel extending from n=9 to n=80. Dissolution of the high lying 
levels has been treated by means of the occupation probability for- 
malism of Hummer & Mihalas (1988), generalised to the non-LTE 
situation by Hubeny, Hummer & Lanz (1994). The calculations in- 
cluded the bound-free and free-free opacities of the H~ ion and 
incorporated a full treatment for the blanketing effects of HI lines 
and the Lyman —a, — /3 and —7 satellite opacities as computed by 
N. Allard (e.g. Allard et al. 2004). For models with T cff >16000K 
we assumed radiative equilibrium while at r e g ^l6000K we in- 
cluded a treatment fo r convective energy tra nsport according to the 
ML2 prescription of iBergeron et al.l dl992D . with a mixing length 
parameter of a=0.6. During the calculation of the model struc- 
ture the hydrogen line broadening was addressed in the following 
manner: the broadening by heavy perturbers (protons and hydro- 
gen atoms) and electrons was treated using Allard' s data (includ- 
ing the quasi-molecular opacity) and an approximate Stark profile 
(Hubeny, Hummer & Lanz 1994) respectively. In the spectral syn- 
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Table 5. Details of the new white dwarf candidate members of NGC 3532 (upper) and NGC 2287 (lower) and the "test" object in the direction of the latter 
cluster (bottom) observed with the VLT. The tabulated spectroscopic signal-to-noise estimates correspond to per resolution element over the range A=4150- 
4300A. Masses a nd cooling times for each star have been estimated using the mixed CO core composition "thick H-layer" evolutionary calculations of the 
Montreal Group (Fontaine et al. 2001). The errors in absolute magnitudes, masses and cooling times shown here are derived by propagating more realistic 
uncertainties in effective temperature and surface gravity of 2.3% and 0.07dex respectively. 



WD 


S/N 


(K) 


logs* 


V 


My 


M 

(M ) 


Tc 

(Myr) 


fj, a cos S 
(mas yr — 1 ) 


(mas yr _1 ) 


Mem? 


Jl 105-585 


120 


13329±^ 


o 1Q +0.06* 


20.33±0.04 


J " L -'°-0.11 


0.72 ±0.04 


370^43 


+4.5±4.1 


-6.9±4.1 


X 


Jl 106-590 


115 


19690::!!! 


o yic+0.03 
S - 45 -0.03 


20.04±0.03 


11 51+ 013 


0.90 ±0.04 


188l» 


-7.9±3.4 


+0.5±3.4 


V 


Jl 106-584 


100 


188331!!® 


o r-9+0.03 
°- oz -0.03 


20.16±0.04 


11 71+0- 13 


0.94 ± 0.04 


2A3+H 


-5.7±3.8 


-1.0±3.6 


V 


J 1107-584 


80 


20923+4°7 




20.19±0.04 


11 69+ 014 


1.00 ±0.04 


210+2g 


-3.9±3.9 


44.6±4.3 




J0645-202 


60 


1 ' ±l 10 _452 


' ■ i " -0.04 


20.97±0.04 


11 26+ - 11 


0.59 ±0.04 


2007^0 


+3.1±4.3 


+9.9±6.4 


X 


J0646-203 


115 




o on+0.06 
S - 8z! -0.06 


20.91±0.03 


u 7 o+0.15 
- LJ "' o -0.1S 


1.12 ±0.04 


™±% 


-4.2±3.4 


-0.3±2.6 


V 


J0644-205 


90 


123767: 233 


7 QH+0-06 


20.52±0.03 


11 ^7+0- 10 


0.59 ±0.04 


334t« 


+6.9±3.3 


-3.4±3.3 


X 



* Formal fit errors - see text for further details. 



Table 6. The co-efficients (zero-point, airmass and colour term) deter- 
mined for Equation 1 , the transformation between instrumental magni- 
tudes and Johnson V. 



K K\ K 2 

27.447±0.020 -0.142±0.016 0.0004±0.0062 



thesis step detailed profiles for the Balmer lines were calculated 
from the Stark broadening tables of Lemke (1997). 



4.2 Effective temperatures and surface gravities 

As in our previous work (e.g. iDobbie et al. 1120091120121) , the com- 
parison between the models and the data is u ndertaken using the 
spectral fitting program XSPEC (Shafer 1991). In the present anal- 
ysis all lines from H-/3 to H-8 are included in the fitting process. 
XSPEC works by folding a model through the instrument response 
before comparing the result to the data by means of a \ 2 — statistic. 
The best fit model representation of the data is found by increment- 
ing free grid parameters in small steps, linearly interpolating be- 
tween points in the grid, until the value of x 2 is minimised. For- 
mal errors in the T c ffS and log gs are calculated by stepping the 
parameter in question away from its optimum value and redeter- 
mining minimum \ 2 urm l the difference between this and the true 
minimum \ 2 c orresponds to \a for a given number of free model 
parameters (e.g. La mpton et al.l 19761) . The results of our fitting pro- 
cedure are given in Table [5] and are shown overplotted on the data 
in Figure [5] 

It should be noted that the formal lcr parameter errors quoted 
here undoubtedly underestimate the true uncertainties. In our sub- 
sequent analysis we have adopted what are considered more realis- 
tic levels of uncertainty of 2.3% and 0.07d ex in effective tempera - 
ture and surface gravity, respectively (e.g. Napiwotzki et aljl99^) . 



WD J1105-5B5 
WD Jl 106-590 
WD J1106-5B4 
WD J1107-5B4 



WD J0645-202 II 

WD J0646-303 1 

WD J0644-305 I 

7 7.5 8 8.5 9 9.5 10 10.5 

m— M 

Figure 6. The derived distance modulii of the white dwarf candidate mem- 
bers of NGC 3532 (upper panel) and NGC 2287 (lower panel) observed 
with the VLT and FORS. The estimated distance modulus of a control ob- 
ject lying in the field of NGC 2287 is also shown (grey typeface). The dis- 
tance modulus of NGC 3532 as estimated by Meynet, Mermilliod & Maeder 
(1993; V-M v =835), Robichon et al. (1999; m-M=8.04) and Kharchenko 
et al. (2005; V-M v =8.6l) and the distance modulus of NGC 2287 as es- 
timated by Harris (1993; V-M V =9A\), Meynet, Mermilliod & Maeder 
(1993; V-M V =9.15), Kharchenko et al. (2005; V-M v =930) and Sharma 
et al. (2006; V-My=930) are overplotted. All new white dwarf candidate 
members have distances consistent with their proposed host clusters, except 
WDJ0645-202 which appears to lie someway behind NGC 2287. 

4.3 Distances 

To place stringent constraints on the distances and to measure the 
proper motions of the seven new spectroscopically confirmed DA 
white dwarfs we have obtained additional observations with the In- 
amori Magellan Areal Camera and Spectrograph (IMACS) and the 
6.5m Magellan Baade telescope. IMACS is a wide field imaging 
system (or multiobject spectrograph) in which light from the Gre- 
gorian secondary of the telescope is fed to one of two cameras, 
either short (f/2.5) or long (f/4.3). The former camera, which was 
used for these observations, contains a mosaic of eight 2048 x 4096 
pixel E2V CCDs. It covers a 27.5' x 27.5' (0.20" pixels) area of sky 
per pointing but suffers from strong vignetting at distances >15' 
from the field center. 
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Figure 7. A vector point diagram of the relative proper motions of the seven spectroscopically observed white dwarfs towards NGC 3532 (left) and NGC 2287 
(right). The proper motions expected to be displayed by members of the two clusters are overplotted (grey open circles). These have been shifted to account 
for the difference between the relative and absolute frames of reference. Objects with proper motions deemed consistent with the cluster populations are 
highlighted (solid symbols). The proper motion of an additional photometric candidate white dwarf member of NGC 3532 is also shown (open triangle, see 
text for further details). 



We acquired 30 s integrations on each object through 
the VbesscII filter during the photometric night of 2010/04/09 
(JD=2455296.25) when the seeing was ~0.7". All targets were cen- 
tered on CCD 2 to negate the effects of chip-to-chip colour sensi- 
tivity variations. The frames were reduced, as per the ESO WFI 
data, using the Cambridge Astron omical Survey Unit CCD reduc- 
tion toolkit Jlrwin & Lewisll200ll) to follow the steps listed in Sec- 
tion [TT] As before, we performed aperture photometry on the re- 
duced images using a circular window with a diameter of 1.5 x 
the full width half maximum of the mean point spread function 
(~0.7"). Several different standard star fields were also observed 
on CCD 2 through both the Vsessett and the /ctio filters during 
the course of this run to allow us to transform our instrumental 
magnitudes onto the Johnson system, 

m = -2.5\og(ADU/t exp ) + K + K X X + K 2 (V - I) (1) 

where ADU is a measure of the total counts from the source, t exp 
the exposure time and X the airmass. The co-efficients and their re- 
spective errors were determined to have the values shown in Table[6] 
and our revised V magnitude estimates for the white dwarfs are 
given in Table[5] For the objects in the field of NGC 2287, this more 
precise photometry is consistent with that based on the ESO WFI 
data. For the NGC 3532 white dwarfs, our refined ma gnitudes are 
in exc ellent agreement with the values determined by Idem et all 
d20 1 lh but typically brighter than our preliminary WFI measure- 
ments by 0.1-0.3 mag. This suggests that our survey of this cluster 
is shallower by 0.1-0.3 mag. at V than was initially concluded. 

We have used our measurements of the effective temper- 
atures and surface gravities of the white dwarfs shown in Ta- 



ble [5] and the model grids o f lBergeronetai1 ( ll995l) . as revised by 
Holberg & Bergeron (2006), to derive absolute magnitudes (My; 
see Table[5]l. Subsequently, we have determined the distance mod- 
ulii of all seven white dwarfs observed with the VLT and FORS, 
neglecting foreground extinction which is low (Av^0.15) along 
the lines of sight to NGC 3532 and NGC 2287. These are plotted 
(solid bars) along with a number of distance estimates available in 
the literature for each of these clusters (dash-dotted vertical lines) 
in Figure |4~3l Based on examination of this figure we consider that 
five out of the six white dwarf candidate cluster members have es- 
timated distances which are consistent with those of their putative 
parent populations (ie. within 2a of the mean). However, we find 
WDJ0645-202 to lie substantially behind NGC 2287 and on this 
basis conclude that it is probably a field object. Not unexpectedly, 
the cooler, redder, "test" white dwarf, WD J0644-205, appears to 
reside in front of this cluster, albeit marginally. 



4.4 Proper motions 

To further explore the cluster membership credentials of these 
white dwarfs, we have exploited the epoch difference of approxi- 
mately a decade between the archival ESO WFI V band datasets 
and our new IMACS V band images to measure their relative 
proper motions. Astrometry is a potentially powerful method of dis- 
crimin anting members of a cluster from the general field population 
(e.g. lHamblvet all 19931) . We have re-examined the images using 
the CASU reduction toolkit to determine the positions, in pixels, 
of bright (but not saturated), unblended, stellar-like objects lying 
within ~2 arcminutes and on the same CCD chip as each candi- 
date. We cross matched these lists of reference star positions using 
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the STARLINK TOPCAT software. Subsequently, we employed 
routines in the STARLINK SLALIB library to construct six co- 
efficient linear transforms between the two images of each candi- 
date, where ^3<j outliers were iteratively clipped from the fits. The 
proper motions, in pixels, were determined by taking the differ- 
ences between the observed and predicted locations of candidates 
in the 2nd epoch imaging. These were then converted into milli- 
arcseconds per year in right ascension and declination using the 
world co-ordinate systems of the datasets and dividing by the time 
baseline between the two observations (~10.12yr for NGC3532 
and 10.35 yr for NGC2287). The uncertainties on these measure- 
ments were estimated from the dispersion observed in the (assumed 
zero) proper motions of stars with similar brightness (jAV|Ss0.5 
mag.) surrounding each white dwarf. The relative proper motion 
vector point diagrams for the objects towards each cluster are 
shown in Figure|7](five-point stars) and the measurements are listed 
in Table[5] The absolute proper motions of each cluster, on th e Hip- 
parcos system dKharchenko et alj|2005l : Ivan Leeuwer]|2009l) , have 
been overplotted (grey open circles) after shifting them by fi a cos 
5, us = +2.8, -1.7 mas yr" 1 and -1.2, -1.8 mas yr" 1 for NGC 3532 
and NGC 2287, respectively, to account for the difference between 
the relative and the absolute frames of reference. These offsets were 
estimated by taking the median of the absolute proper motions of 
stars at the faint end of the UCAC3 catalogue which reside within 
the survey areas. We note that there is a substantial overlap with 
our reference stars. 

Of the three white dwarfs in the direction of NGC 2287 that 
were observed with the VLT and FORS2, we find only WD J0646- 
203 (solid star in right hand panel) has a proper motion which we 
deem to be consistent with this cluster (ie. within 2a). While three 
objects (solid stars in left hand panel) towards NGC 3532 have 
proper motions within la of the cluster mean value, WD Jl 105-585 
(left most open star) lies more than 3 it away and on these grounds 
seems unlikely to be associated with NGC 3532. The proper mo- 
tion of NGC 3532- WDC Jl 105-5842 (open triangle), an additional 
object which we were able to include in our astrometric study due 
to its spatial proximity to WD Jl 105-585, is found to be discrepant 
from that of the cluster at >4a (p a cos <5, /j<5=-20.9±3.0,+2.8±5.4 
mas yr" 1 ). This object has similar colours to WD Jl 106-590 and 
WD Jl 106-584, B-V^0.l and despite lying close enough to the 
theoretical tracks in the B-V, V CMD to be considered a photo- 
metric cluster member, was too bright (V~19.7) to have featured 
in our spectroscopic follow-up program. 



Table 7. Inferred progenitor masses of all the known white dwarf mem- 
bers of NGC 2287 for the most likely cluster age (column 2) and plau- 
sible limiting ages (columns 3 and 4). The cooling times of all white 
dwarfs were initially derived using CO core tracks. The cooling time 
of WD J0646-203 was also evaluated using ONe models (bottom row). 
The estimated stellar lifetimes have been t ransformed to mass es using 
the solar metallicity evolutionary models of Girardi et al. 1 2000). 



5 NEW SOUTHERN OPEN CLUSTER WHITE DWARFS 

The four objects identified in Section|4]as having properties consis- 
tent with cluster membership represent the faintest white dwarfs in 
NGC 3532 and NGC 2287 unear thed to date. Once a gain we have 
used the evolutionary models of lFontaine et al.l d200ll) to determine 
their masses and cooling times. For completeness we have also 
derived the masses and cooling times for the three white dwarfs 
in our spectroscopic sample which we have already rejected as 
probable field objects. So that this work is consistent with our 
previous efforts and those of several other recent studies in this 
area ( e.g ICasewell et al.l|2009l; iDobbie et alj|2009l : Iwilliams et all 
2009; Kalir ai et al J 120081) we have (initially) adopted the calcula- 
tions which include a mixed CO core and thick H surface layer 
structure. 



ID 




M init (M ) 




Adopted cluster age 


243Myr 


203Myr 


283Myr 


NGC 2287-2 


4 45+0-20 


4 99+<>-28 
^■' v -0.22 


4.08+.^ 


NGC 2287-5 


4 57+0.24 


5 16+ 037 


4.1712:1! 


WDJ0646-203 


g 41 + 1.54 


9- 73 -2.58 


c 99 +0.68 

-0.39 


WDJ0646-203 




8.06+J; 1 * 


5 02+ - 51 

J - u -0.32 



5.1 Is WD J0646-203 in NGC 2287 a ONe core white dwarf? 

Based on our measurements of effective temperature and surface 
gravity we have determined the mass of WDJ0646-203 to be 
M=1.12±O.O4M (Table [5}. While the astrometric characteris- 
tics, the distance modulus and the projected location of this white 
dwarf only ~0.3 7° from the cluster center (ie. wit hin the projected 
tidal radius e.g. ICoxlll954t iPiskunov et alj 2008) are firmly sug- 
gestive of an association with NGC 2287, we have nonetheless ex- 
plored the likelihood that it is merely a field star. At T ef r>25000K, 
WD J0646-203 is hot enough to radiate significantly in the extreme- 
ultraviolet (EUV) and, had it been located within the Local Bub- 
ble, would likely have been detected i n the EUV al l -sky s urveys 
dPve et al.1 1 19951 : iBowver et all 1 19961) . IVennes et alj i 19971) have 
computed the space de nsity of EUV detect ed white dwarfs to be 
n^l.6-2.2xl0" 5 pc" 3 . iLiebert et all j2005h have estimated that 
only ~l/5th of white dwarfs in the Palomar-Green survey have 
M>0.8M Q while lVennes et alj d 19971) find that 10 out of their sam- 
ple of 90 have M<;1.1M . We estimate that as a result of flagging 
as potential members those objects with distance modulii consistent 
with the cluster mean at the 2a level, our survey has probed a vol- 
ume of ~9500pc 3 . Thus, with the above space densities we should 
have expected to unearth ~0.04 hot, massive field white dwarfs in 
our study of NGC 2287. Assuming a poisson distribution, we de- 
termine that the probability of us detecting one field object with 
the characteristics of WD J0646-203 is rather low, P^0.04, and we 
are led to conclude that this white dwarf is a bona-fide member of 
NGC 2287. 

Some theoretical studies suggest that the surface H-layers 
of the most massive white dwarfs may be of significantly lower 
mass than t he 10~ 4 Mq adopted in the "standard" evolutionary 
models (e.g. iGarcia-Berro et al1ll997l) . Consequently, we have re- 
evaluated the parameters of WD J0646-203 using the thin H-layer 
models of the Montreal group which allow for a surface H mass 
of only 1O~ 1O M . In this case we obtain a slightly lower mass 
and marginally shorter cooling time but the differences between 
the two determinations are less than O.O2M and 3Myr, respec- 
tively. Since we know the effective temperature of WD J0646-203 
from our spectroscopy, to perform something of a sanity check on 
these mass determinations, we have also estimated it by assuming, 
a priori, that this white dwarf lies at the distance of the cluster, cal- 
culating its radius and applying a theoretical mass-radius relation. 
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Table 8. Inferred progenitor masses of all the known white dwarf mem- 
bers of NGC 3532 for the most likely cluster age (column 2) and plau- 
sible limiting ages (columns 3 and 4). The white dwarf cooling times 
were derived using CO core tracks and the stellar lifetimes have been 
tr ansformed to masses using the solar metallicity evolutionary models 
of lGirardi etail feOOOh . 



ID 




M init (M ) 




Adopted cluster age: 


300Myr 


275Myr 


325Myr 


NGC 3532-1 


, „t+0.09 
J - < "-0.07 


■j qq+0-10 
-'•"-0.09 


3 69+ 08 
J - Dy -0.06 


NGC 3532-5 


o 71 +0.06 

■'•"-o.os 


oc+0.07 
J -° -0.06 


* cn+0.06 

J - J '-0.04 


NGC 3532-9 


, <; 7 + 0.02 

J - J '-o.oi 


3 69+ 02 


3 47+0.01 

J - H '-o.oi 


NGC 3532-10 


45 „+0.35 


4 Q2+ - 43 
^ z -0.31 


4 30+ 28 

^■ JU -0.19 


WD Jl 106-590 


c ]4 +0.69 
J - 1 -0.39 


570 +i.i9 


4 76+ 047 

4 -'°-0.32 


WD Jl 106-584 


6.92+« 2 8 




S Sfi +1 ' 75 
3 - St) -0.75 


WD Jl 107-584 






5 ll+ - 79 
J - 11 -0.42 



If the assumption of cluster membership is safe, the systematic un- 
certainty in the mass estimate, which relates to our spectroscopic 
analysis method and our choice of at mosphere code for generat- 
ing the synthetic H line profiles (e.g. Lie bert et al.ll2005l) . should 
be reduced here. Using the effective temperature shown in Table[5] 
the arithmetic mean of the four distance modulii listed in the cap- 
tion of Figure H31 (my -Mv=9. 29+0.09) and the mass-radius re- 
lation from the thick and thin H-layer CO core evolutionary mod- 
els we determine WDJ0646-203 to have Af=l.O9±O.O3M and 
Af=l.O7±O.O3M , respectively. 

These values argue that WD J0646-203 is potentially the most 
massive white dwarf yet identified within an open cluster. For 
example, LB 1497, the sole white dwarf located within the pro- 
jected tidal radius of the Pleiades cluster, has M= l .Q5±O.Q2M , 
based on a weighted mean of the estimates in iBergeron et alj 
dl995h . iDobbie et al.1 d2006h and lVennes & Kawkai l2008h . which 
were derived using thick H-layer CO core evolutionary models. The 
masses of the four most massive white dwarfs recently unearthed 
in the slightly metal poor cluster NGC 2168, which were also de- 
rived using thick H-layer CO core tracks, appe ar to be clumped 
around Af =1.O1-1.O2M dWilliams et ai]|2009l) . Indeed, our CO 
core based mass determinations for WD J0646-203 sit on or above 
recent theoretical predictions for the minimum mass of ONe core 
white dwarfs (e.g. M~1.O5-1.1M , ISiessI 120071 ; iGil-Pons et all 
l2003h . The most likely progenitor mass we infer of A/j n it~6.3- 
6.4M (thin or thick H-layer white dwarf models) is also consistent 
with the initial mass range where some solar metallicity stellar evo- 
lutionary calculations, which include overshooting from the con- 
vect ive regions of the star, predict central C igni tion (Afi n it>6M 
e.g. lEldridge&Tou3l2004t iBertelli et all 120091) . While the igni- 
tion of C in the core is anticipated to occur only at larger initial 
masses (by ~1-2M ) in stellar models that neglect to treat this 
phenomenon, a number of distinct observational based results pro- 
vide compelling arguments in favour of significant convective over - 
shooting (e.g. (woo et al.l 120031 : IClaretJ 120071 : ISalaris et al.l l2009h . 
Moreover, despite our inferred progenitor mass ( A/i n jt~5.2M ) 
for a cluster age towards the upper end of the plausible range for 



NGC 2287 (see Table|7j>, probably being too low for a star of this 
composition to have experienced central C burning, in this eventu- 
ality, WD J0646-203 would appear as something of an outlier in a 
semi-emprical IFMR, sitting iJO. 1 M above all other white dwarfs 
from solar metalicity populations located in this pa rt of the initial 
mass domain. For exampl e, the well studied Sirius B dBarstow et al.l 
l2005t iLiebert et al.ll2005l) and the NGC 3532 and NGC 2516 ob- 
jects studied by us here using a similar spectroscopic set-up and 
the same models and analysis techniques cluster around Af~0.95- 
l.OM . Of course, this variation could be statistical in nature but 
earlier we noted that our estimates of mass obtained using two dif- 
ferent approaches are in good agreement. 

Bearing in mind that objects of ONe composition are ex- 
pected to adhere to a slightly different mass-radius relation, we 
have re-evaluated the parameters of WDJ0646-20 3 once again, 
this tim e using the ONe core evolutionary tracks of lAlthaus et al.1 
J2007I) . Based on our spectroscopic effective temperature and sur- 
face gravity measurements we estimate a mass and a cooling time 
of M=1.O9±O.O4M and r coo i=164^^Myr, respectively. Follow- 
ing our other approach, where we assume WD J0646-203 lies at the 
distance of NGC 2287, we determine M=l .05±0.03M Q . These re- 
vised parameters for our object remain in accord with modern the- 
oretical estimates of the minimum mass of an ONe core UMWD 
at approximately solar metallicity dSiessll2007h and, for the most 
likely cluster age, lead to an inferred progenitor mass consistent 
with the initial mass range where central C ignition is anticpated 
(see Table [7}. Clearly it not possible for us to definitively deter- 
mine the composition of the core of WD J0646-203 but on the basis 
of our new observations and guided by modern theoretical calcula- 
tions, we conclude that there is a distinct possibility that it consists 
of O and Ne, having evolved from a single, heavy-weight, interme- 
diate mass progenitor star. 

5.2 The NGC 3532 white dwarfs 

In our previous work on NGC 3532, we revisited six white dwarfs 
in the direction o f the cluster and confirmed four as members 
dPobbie et al.ll2009h . Two objects, NGC 3532-6 and NGC 3532-8, 
were rejected by us on the basis that our revised estimates of their 
distances are incompatable with that of the cluster. We have now 
found three further white dwarfs members of NGC 3532 taking the 
total number to seven and making this is one of the richest hauls 
of cluster degenerates where good quality optical spectroscopy has 
been secured. This population is important since NGC3532 has a 
turn-off mass of Mi n i t «3.5M and spans the age gap between the 
younger NGC 2168 and more mature Praesepe and Hyades clus- 
ters. 

Despite our new study probing to fainter magnitudes than the 
previous photographic survey of the cluster, we have not identi- 
fied any compelling candidate ultra-massive ONe core white dwarf 
members of NGC 3532 (ie. M^1.O5M and A/ ini t^6M ). In- 
stead, the three new fainter degenerate members have masses that 
are similar to that of the previously known faintest white dwarf 
in the cluster, NGC 3532-10. Intruiguingly, despite their compa- 
rable masses, we find that the cooling times of these four most 
massive white dwarf members span a range of ~90-100Myr, ir- 
respective of our choice of thick or thin H-layer evolutionary mod- 
els. For example, for WD Jl 106-584 we estimate r coo i~240Myr 
while for NGC 3532-10 we determined r coo i~150Myr. We sug- 
gested in our previous study, where we analysed an ensemble of 
~50 white dwarfs from ~10 mainly cluster populations, that the 
IFMR is less steep at Afi n i t «4M than immediately below this 
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Figure 8. The location of the seven known NGC 3532 white dwarfs in initial mass-final mass space for the most likely cluster a ge (300Myr, black triangles) 
and approriat e limiting ages (275Myr, grey circles and 325Myr, grey squares). A theoretical IFMR from Marigo & Girardi ( 2007) (dark-grey dashed line) and 
a 1TPR from Wagenhuber & Groenewegen ( 1998) (light-grey dot-dashed line) are overplotted. 



initial mass regime. We noted that the distribution of these white 
dwarfs mirrored the form of the theoretical relationship between 
the initial mass of a star and the mass of the CO core at the time 
of the first therma l pulse (1TPR) as it evolves onto th e asym ptotic 
giant branch (e.g. iBecker & Ibenlll979t) . ISalaris et all d2009t) inde- 
pendently reached a similar conclusion regarding the slope of the 
(semi-) empirical relation by comparing the heterogeneous open 
cluster data to an IFMR predicted by the BaSTI stellar evolution- 
ary models of lPietrinferni et al.ld2004l) . 

Nevertheless, given the rather heterogeneous nature of that 
ensemble of white dwarfs, it is conceivable that cluster to cluster 
age uncertainties could have conspired in such a way to mislead 
us. For example, if the ages of the youngest populations which 
dominate the upper end of the relation are systematically under- 
estimated with respect to the more mature populations which dom- 
inate the lower end then the semi-empirical IFMR would appear 
less steep at higher initial masses. However, independent of both 
the age we adopt for NGC 3532, within plausible limits, and our 
choice of white dwarf evolutionary model grid (thick or thin sur- 
face H-layers), the location in initial mass-final mass space of these 
seven degenerates (Figure [8] and see also Table [8), which span a 
broad range in initial mass yet are vulnerable to only a single clus- 
ter age unc ertainty, are also consistent with the t rend outlined by the 
1TPR (e.g. Wagenhub er & Groene wegen 1998) and the theoretical 
IFMR (e.g. Marig o & Girardill2007l) . This lends some further sup- 
port to both our earlier conclusion and that of ISalaris et al.l (2009) 
regarding the form of the IFMR. 

Population synthesis calculations which adopt an IFMR with 



a lower gradient above Mi n it«4MQ, than in the immediately 
lower initial mass regime, better reproduce the shape of the field 
white dwarf mass distribution, in particular the secondary peak at 
M^0.8-1.0Mp), than those which e mploy a relation with a sim- 
ple linear form dFerrario et al.ll2005l) . It has been proposed previ- 
ously that <;80% of white d warfs with M>Q,8 Mq might originate 
from binary mergers (e.g. lLiebertetal.1120051) but such a change 
in the gradient of the IFMR mitigates the need to invoke such an 
evolutionary channel to explain their existence. That stellar evo- 
lution might be able to readily produce large CO cores without 
the need for mergers is potentially important in the context of 
Type la supernovae which are believed to result from the com- 
plete thermonuclear disr uption of massive CO white dwarfs (e.g 
lYungelson & Liviol200oT) . For example, in a leading, current, theo- 
retical formation channel, the sub-Chandrasekar double-detonation 
model, the more massive a CO core progenitor, the less material is 
required to be accreted to trigger the explosion. We note that the 
lifetimes of those stars to which the less steep section of the IFMR 
is applicable (ie. Mt n it^4M0) are comparable to the timescales of 
the prompt componen t of the Type la delay time distribution (e.g. 
see lRuiteret alfeoill) . 



6 SUMMARY 

We have performed a CCD photometric survey, complimented by 
optical spectroscopic and astrometric follow-up studies, of ~0.3 
square degrees of sky towards each of the young/intermediate age 
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open clusters NGC3532 and NGC2287, in a bid to unearth in- 
trinsically faint white dwarf members. This work has led to the 
identification of four new white dwarfs which are probably clus- 
ter members, three in NGC 3532 and one in NGC2287. Spectro- 
scopic analysis of these objects has revealed WD J0646-203 to be 
potentially the most massive white dwarf member of an open clus- 
ter unearthed so far, with there being a distinct possibility of it hav- 
ing a core composed of O and Ne (Af =1.02-1. 16Mq). Our study 
has also shown that, despite a range of ~90Myr in age, the four 
most massive degenerates in NGC 3532 are rather similar in mass, 
consistent with our earlier suggestions that the IFMR could be less 
steep at Mi a it^4MQ than at initial masses immediately below this 
(Figure [8}. The results of our analysis of new and improved spec- 
troscopy of the four known NGC 25 16 white dwarfs are also in 
agreement with our prior conclusion. 

However, further studies of young/intermediate age, rich, open 
clusters remain necessary to better determine the relative frequen- 
cies of massive and ultra-massive degenerates and to pin down the 
maximum mass of a white dwarf that can form via canonical stel- 
lar evolution. It appears probable that additional constraints on the 
form of the IFMR will come from NGC 3532 since a more exten- 
sive CCD survey of the cluster has re cently unearthed m any more 
faint candidate white dwarf members l iClem et alll201 lb . Low res- 
olution optical spectroscopy with good signal-to-noise is now re- 
quired to characterise these objects. Substantial amounts of time 
on the current generation of large telescopes would be required to 
obtain spectroscopy of the faintest white dwarf members of attrac- 
tive, even richer but more distant clusters. Alternatively, it might be 
possible to address this type of investigation photometrically since 
simulations suggest that in the U,U-V colour-magnitude diagram 
of a very rich cluster with r~200-300Myr (e.g. NGC 6705), a pu- 
tative ultra-massive white dwarf population would define a rather 
prominent blue hook at the bottom end of the cooling sequence and 
give rise to a secondary peak in the U band luminosity function. 
Detection of this feature would need photometry of good quality 
down to Vf»25 so even this approach could prove rather challeng- 
ing given the likely levels of scattered light from the large number 
of intrinsically bright stars in rich clusters of this age. 
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